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ABSTR.\CT 


An application of the NASA unstructured grid software 
system TetrUSS is presented for the prediction of 
aerodynamic drag on a transport configuration. The paper 
briefly describes the underlying methodology and 
summarizes the results obtained on the DLR-F4 transport 
configuration recently presented in the first ATAA CFD 
Drag Prediction Workshop. TetrUSS is a suite of loosely 
coupled unstructured grid CFD codes developed at the 
NASA Langley Research Center. The meshing approach is 
based on the advancing-front and the advancing-layers 
procedures. The flow solver employs a cell-centered, finite 
volume scheme for solving the Reynolds Averaged Navier- 
Stokes equations on tetrahedral grids. For the present 
computations, flow in the viscous sublayer has been 
modeled with an analytical wall function. The emphasis of 
the paper is placed on the practicality of the methodology 
for accurately predicting aerodynamic drag data. 


INTRODUCTION 


Computational Fluid Dynamics (CFD) has matured to 
the point where quality aerodynamic information can be 
routinely produced for many engineering applications. The 
level of success in implementing CFD for solving realistic 
problems depends on several factors. Obviously, the 
accuracy of the underlying methodology directly affects the 
viability of solutions and is of primary importance. 
Moreover, the robustness of the CFD method and the 
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convenience with which it can be employed determine its 
productivity for solving realistic problems. The timeliness 
of generating numerical data is often as significant as 
accuracy of the solutions, both contributing to the 
usefulness of the final result. The challenge is especially 
increased in a design environment where a large amount of 
data must be generated and analyzed in a short period of 
time. Therefore, a CFD capability must not only meet the 
quality requirement of the predicted infonnation but also 
the practicality and automation with which the data are 
mass-produced, for example, for generating extensive drag 
polar and drag rise data. 

Although recent years have witnessed many 
breakthrough advances in the related computational fields, 
the suitability of CFD methods have not been fully 
evaluated for many real-world industrial applications. One 
of the critical challenges of CFD for predicting commercial 
aircraft perfonnance is the ability of the method to produce 
reliable drag data, Aerodynamic drag prediction and 
reduction remains an important (and often fonnidable) task 
for the commercial aircraft design to reduce operational 
costs. While progress in various aspects of CFD are 
achieved and reported on a regular basis, less attention has 
been focused on the subject of drag prediction and 
validating the current CFD methodologies for such 
applications. 

In that regard, a workshop was recently organized by 
the AIAA Applied Aerodynamics Technical Committee to 
address the issues related to CFD drag prediction for 
commercial aircraft design and to evaluate the existing 
computational capabilities for this purpose. The objective 
of the workshop was, in particular, to assess the state-of- 
the-art computational methods as practical aerodynarnic 
tools for aircraft force and moment computation with 
emphasis placed on prediction of drag. The workshop was 
held in conjunction with the 19*^ AIAA Applied 
Aerodynamics Conference in Anaheim, CA, 9-10 June 
2001 and was attended by several international participants 
who presented their results on a common wing/body 
configuration. The workshop provided an impartial forum 
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for evaluating the effectiveness of the existing CFD codes 
as applied to the subject of drag prediction. 

The geometry selected for the workshop activity was a 
subsonic wing,T)ody transport configuration referred to as 
DLR-F4 model. To maintain as much uniformity among 
various computations as possible, several baseline multi- 
block structured, overset, and unstructured grids (similar in 
resolution) were generated by the workshop organizing 
committee and provided to the participants to perform their 
baseline (required) computations. Participants were also 
encouraged to develop their own grids and apply their best 
practices to perform additional (optional) cases. The results 
generated by the participants are included in the workshop 
proceedings^ A summary of the presented data is also 
given in a companion paper by the workshop organizing 
committee^. 

This paper serves as an extended documentation for the 
baseline unstructured grids provided to the workshop 
participants and also for one of the participating CFD 
activities. The paper briefly describes the underlying 
unstructured grid methodology and presents the results 
obtained for the DLR-F4 configuration for a range of Mach 
numbers and angles of attack. 


MODEL AND FLOW DESCRIPTION 


As stated in the workshop summary paper', several 
requirements were considered by the selection committee 
for choosing an appropriate configuration for the workshop 



Figure I. DLR-F4 wing/body geometry (Ref. 3). 


computational activity. Among the selection criteria, the 
following were given special consideration: 1) relevancy to 
the workshop’s interest in the transport aircraft design for 
transonic flight, 2) availability of reliable experimental 
data, 3) relative geometric simplicity for ease of performing 
computations by all workshop participants, and 4) non- 
classified, non-proprietary requirements. After careful 
examination of several options, the DLR-F4 wing/body 
configuration was selected, as it met all of the above 
conditions. The model has been extensively tested in three 
different European wind tunnels: ONERA-S2MA, NLR- 
HST, and DRA -8ftX8ft DRA Bedford\ It has also been 
utilized for CFD code validation in Europe"*. The lest results 
are freely available and, along with the geometry 
description, are documented in detail in Ref. 3. 

The subject geometry consists of a fuselage 1 192 mm 
long and a swept-back wing of 585.7 mm semi-span. The 
wing aspect ratio is 9.5 with a leading-edge sweep angle of 
27.1 degrees (see Fig. 1). The wing tip is defined by half- 
circular cross-sections with local diameters equal to the 
corresponding wing thickness values at the tip. The 
wing/body junction is sharp, identified by only the 
intersection curv'es (no fillet defined). The body and the 
wing are numerically specified in Ref 3 by a number of x- 
and y-constant cross sections, respectively. Although 
geometrically simple, the model resembles a modern 
transport-type aircraft configuration, which along with 
extensive experimental data, provides a good candidate data 
base for CFD code validation. 

Using the given numerical sections, the workshop 
organizing committee constructed a Computer Aided 
Design (CAD) surface definition by linear lofting between 
defining sections. The CAD definition was used for 
generating the baseline CFD grids and also for distribution 
among the workshop participants. The CAD definition, 
given in standard Initial Graphics Exchange Specifications 
(IGES) format, contains 45 surfaces and 4 trim curves 
shown in Figure 2. For computational economy, a semi- 
span model of the geometry is constructed. The wing 



Figure 2. Surface IGES definition for DLR-F4 
wtng/body geometry. 
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surface definition incorporates an aeroelastic deformation 
under load based on a theoretical estimate for the design 
condition reported in Ref. 3. 

The DLR-F4 model is designed for flying at subsonic 
through transonic speeds. The dominant cruise flow 
condition, as observed from the wind tunnel experiments, is 
primarily attached transonic flow over the wing with 
pockets of supersonic llow' on the upper surface tenninated 
by weak shock waves. At the design condition, a small 
separation of flow is observed at the kink region of the 
wing trailing edge. All experimental data given in Ref. 3 
are corrected to "Tree air” conditions. 


tetrahedral cells grow in the computational field from a 
triangular surface mesh (initial front). The advancing 
process continues until the entire domain is filled with 
contiguous tetrahedral cells. The advancing- front technique 
inserts individual cells in the inviscid portion of the flow 
field in an irregular fashion. The lack of an apparent order 
in which the grid is generated contributes to the flexibility 
of the method. On the contrary, the advancing-Iayers 
method generates thin layers of tetrahedral cells packed in 
the boundary-layer region in a more orderly manner. Each 
advancing layer contains a stack of three tetrahedral cells. 
The systematic way of generating grids by the advancing- 
Iayers method is favored for the ease of generation and 
better quality of thin tetrahedral cells in the boundary layer. 


COMPUTATIONAL METHOD 


The CFD technique employed for the present 
computations is based on a tetrahedral unstructured grid 
methodology. Unstructured grids have gained increasing 
popularity in the CFD community mainly due to two salient 
features: 1) the inherent flexibility of the grid elements for 
discretizing complex domains and 2) the convenience with 
which the generated grids can be post-processed and 
modified. While the fonner reduces the computational 
cycle time by nearly an order of magnitude for a typical 
complex configuration, the latter facilitates the 
implementation of automatic grid adaptation to the flow 
and/or geometric features. 

A coordinated activity has been underway at the NASA 
Langley Research over the past decade to develop an 
integrated system of unstructured grid codes. The primary 
objective of this team effort was to bring the state-of-the-art 
in CFD to a higher level of usability in the analysis and 
design environment. The outcome has been a system of 
user-friendly software referred to as TetrUSS^ The system 
consists of a tetrahedral grid generation package (GridTool 
and VGRTDns), an Euler and Navier-Stokes solver 
(USM3Dns), a post-processing analysis code (ViGPLOT), 
and several other tools and extensions of the codes for 
solving specific problems. An adaptive grid capability has 
also been devised recently and is under further development 
at present time.^ Although these codes are often used in 
combination, the TetrUSS system is modular allowing 
external unstructured grid codes to be used within the 
system synergistically. 

A brief description of the TetrUSS components is 
presented below as applied to the DLR-F4 case for 
completeness. Further information about the underlying 
methodology can be found in the respective references. 


Grid Generation 

The unstructured grid code VGRIDns generates single- 
block, tetrahedral grids for both Euler and Navier-Stokes 
computations. The grid generation method is based on the 
Advancing-front^ and the Advancing-Iayers^ techniques. 
Both techniques resemble marching procedures by which 


The generation of a surface mesh is accomplished in 
two main steps: 1) geometry set-up and preparation of grid 
parameters, both using the graphical grid utility code 
GridTool^ and 2) surface triangulation with VGRTDns. The 
geometry of interest is first defined in terms of a set of bi- 
linear parametric and/or Non-Uniform Rational B-Spline 
(NURBS) patches using the given CAD surface definition. 
The surface patches are constructed with GridTool 
interactively. The union of all these contiguous patches 
(including those defining the outer boundaries) forms a 
"water-tighC surface enclosing the entire computational 
domain. Figure 3 shows a surface patch representation of 
the DLR-F4 model for generating unstructured grids with 
VGRIDns. In this example, the model, along with a 
rectangular outer boundary (not shown), is defined in terms 
of 106 bi-linear parametric and NURBS patches. 

The distribution of the unstructured grid points on the 
surface and in the field is accomplished with VGRIDns by 
means of a set of user-prescribed source elements. These 
sources are also defined and placed in the domain using 
GridTool. The source elements, along with a “transparent” 
background grid^^ (on which the source information is 
distributed by solving a Poisson equation), provide the 
required data for generating the final unstructured grid. 





Figure 3. Surface patch representation of the DLR-F4 
model. 
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Figure 4. Arrangement of source elements on DLR- 
F4 for grid clustering. 


Figure 4 illustrates several “linear” source elements placed 
in the DLR-F4 fuselage and wing for controlling the grid 
clustering on and around the geometry. 

The process of preparing surface patches and grid 
distribution parameters constitutes 50-90% of the total grid 
generation time depending on the complexity of the 
geometry and the "cleanliness'’ of the CAD surface 
definition. In any grid generation method, a substantial 
amount of time and effort is often spent to repair faulty or 
incomplete CAD surface definitions before the actual 
process of grid generation starts. For the present problem, 
the surface geometry was well defined, allowing the surface 
patches to be constructed easily in approximately one hour. 
The complete first grid, from CAD to the final post- 
processed volume grid, was generated in one day on a 
Silicon Graphics (SGI) workstation. Additional iterations 
were required to modify source parameters to fit the final 
grid with the desired characteristics (such as size, local 
resolution, stretching, etc.) prescribed by the committee for 
the workshop computations. 

Once the surface patches and sources are in place, the 
surface mesh is generated automatically with VGRIDns. In 
VGRIDns, individual surface patches are triangulated using 
the two-dimensional (2D) versions of the advancing-front 
and advancing-layers methods. Additionally, a recent 
modification of the code allows triangulation of the 
NURBS surfaces in the 3D space. The new capability 
provides for better quality surface grids, guarantees the grid 
nodes to stay on the CAD surface, and accepts surface 
patches of arbitrary shapes and number of edges (a 
substantial simplification of the patching process). The 
capability of direct surface triangulation on NURBS 
surfaces is still under further development and is not 
available in the current (released) version of VGRIDns. 



Figure 5. Surface triangulation on the DLR-F4 model 
and the plane of symmetry . 


A collection of all triangulated surface patches makes 
up the final surface mesh that also acts as an initial front for 
generating the tetrahedral volume grid by the advancing- 
front/layers methods. Figure 5 shows the surface 
triangulation on the DLR-F4 model as well as the grid on 
the symmetry plane. 

A salient feature of VGRIDns is a multi-directional, 
anisotropic grid clustering capability, by which grid points 
can be distributed nonuni formly in different directions, 
producing stretched grid elemen^^^^ result is a reduction 
of total cell count by a factor of three without losing the 
grid resolution in the essential directions. For the present 
problem, grid stretching is mainly applied at the wing 
leading and trailing edges in the spanwise direction while 
maintaining high concentration of grid nodes chordwise 
(see Fig. 6). 

As mentioned earlier, the generation of tetrahedral 
volume grid is achieved with the 3D versions of the 
marching techniques, similar to those utilized for surface 
triangulation. After the process of volume grid generation 
with VGRIDns is over, some voids may remain in the 
inviscid portion of the grid which cannot be closed with the 
grid generator itself. Furthermore, the generated grid 
usually contains isolated pockets of distorted cells, which 
may be unacceptable to some solvers. The inherent 
irregularity of unstructured grids not only provides ease of 
generation but also offers flexibility for alteration and post- 
processing of the generated grids. Any segment of a 
tetrahedral grid can be removed and locally re-meshed 
without disturbing the rest of the grid. This operation is 
perfonned with a computer program referred to as 
POSTGRID as the last stage of the generation process. 
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Figure 6. Details of surface grid on the DLR-F4 model: (a) wing upper surface showing stretched grid in 
spanwise direction (b) variable stretching in the wing tip and leading edge region, (c) ‘‘full- 
viscous” grid with 35 viscous layers, and (d) “wall-function” grid with 12 viscous layers. 


Flow Solution 

USM3Dns is a tetrahedral cell-centered, finite volume 
Euler and Navier-Stokes solver^'. The inviscid fiux 
quantities are computed across the cell faces using the Roe's 
Hux-difference splitting scheme, and the spatial 
discretization is accomplished by a novel analytical 
reconstruction process. The solution is advanced in time to 
steady state using an implicit backward-Euler time-stepping 
scheme. Flow turbulence is modeled by the Spalart- 
Allmaras one-equation model, which is optionally coupled 
with a wall function to reduce solution stiffness and the 
number of cells in the sublayer of the boundary layer. All 
computations presented in this paper are performed using 


the wall-function feature of USM3Dns in a fully turbulent 
mode. 

USM3Dns runs on massively parallel computers and 
clusters of personal computers (PC's). A single processor 
version is available for vector processors, such as the Cray 
super-computers, with multi-tasking, and single 
workstations. The parallel version is the version of 
choice because it enables rapid turnaround of large 
problems. The code requires 175 eight-bit words of 
memory per tetrahedron. It runs with a speed of 34 
psec/cell/cycle/processor on a Cray C90 and 230 
psec/cell/cycle/processor on the SGI Origin 2000 parallel 
computer. 

USM3Dns supports standard boundary conditions such 
as the flow tangency, no-slip solid surface, characteristic 
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inflow/outflow {for subsonic flows), and freestream- 
inflow/extrapolation-outflow (for supersonic flows). In 
addition, some special boundary conditions including wall 
functions, wake flow, jet engine, and propeller are available 
in the code, and are reviewed in Ref 5. 


RESULTS 


Tetrahedral Grids for Workshop 

For the required workshop CFD computations, three 
unstructured grids were generated on the DLR-F4 wing- 
body configuration and distributed among the workshop 
participants. Table I contains the specifications of these 


grids. While the basic characteristics of the grids (e.g., 
“effective’’ resolutions) are comparable around the 
geometry, each grid has been generated for a specific 
solution strategy. Because tetrahedral unstructured grids 
contain 5-6 times more cells than nodes, a grid of sufficient 
density for a cell-centered solver is usually considered too 
coarse for a solution technique utilizing a node-based 
scheme. In addition, solvers employing analytical wall 
functions to approximate flow in the viscous sublayer 
require fewer grid elements across the boundary layer 
compared to those solving the Navier-Stokes equations to 
the surface. 

The distributions of grid spacing normal to the wall 
were prescribed uniformly over the entire geometry. The 
nominal nonnal spacings were determined at the mid-span. 


Grid Characteristics 

Node-based Grid 
(Full Viscous^ 

Cell-based Grid 

Full Viscous 

Wall- function 

Boundary points 

48,339 

23,290 

25,175 

Surface triangles 

96,674 

46,576 

50,346 

Triangles on no-slip surfaces 

72,902 

30,037 

38,571 

Total grid points 
Points in viscous layers 

1,647,810 

470,427 

414,347 

1,129,427 

389,753 

238,301 

Total tetrahedral cells 


2,743,386 


Cells in advancing (viscous) layers 


2,208,260 


Vlaximum number of viscous layers 

35 

35 

12 

Number of complete viscous layers 


24 

7 

Grid points across wing T.E. 


5 

5 

Chordwise grid spacing at L.E. 


-- 0.450 mm 

~ 0.450 mm 

Chordwise grid spacing at T.E. 

-- 0.500 mm 

- 0.800 mm 

- 0.800 mm 

Maximum spanwise spacing at L.E 

~ 2.500 mm 

~ 6.000 mm 

- 6.000 mm 

Maximum spanwise spacing at T.E 

~ 3.500 mm 

~ 3.500 mm 

- 3.500 mm 

Grid spacing on the fuselage 

~ 10.000 mm 

- 10.000 mm 

~ 10.00 mm 

Grid spacing at the outer boundary 

- 3000.000 mm 

3000.000 mm 

~ 3000.000 mm 

First “viscous” spacing off the wall 

0.001 mm 

0.003 mm 

0.0549 mm 

Rate of geometric stretching 
(viscous layers) 

1.2 

~ 1.2 

Eq.(l) 

Outer boundary box size 

50 mean chords 

50 mean chords 

50 mean chords 


Table 1. Grid specifications for three different unstructured grids generated for the DLR-F4 Drag 
Prediction Workshop. 
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mid-chord location on the wing using flat-plate turbulent 
boundary layer theory. For the ‘‘wall-function” grid, an 
initial spacing for of 50 was assumed, whereas those of 
the “full viscous” Navier-Stokes grids were set to provide a 
of 1. For the full Navier-Stokes grids, a geometric 
stretching factor of 1,2 was applied. For the wall-function 
grid, the clustering of points normal to the surface was 
computed according the stretching function 

S ,- 5,1 1+^(1 ( 1 ) 

where 5„ is the normal spacing of the layer, 8, is the first 
layer spacing prescribed by the user, and the factors a and b 
are constants that determine the rate of stretching, also 
prescribed by the user. For the current wall-function grid, 
the parameters 5,, a, and b have been set to 0.0549 mm, 
0.37 and 0.07, respectively, which have produced 
approximately 18 tetrahedral cell layers (6 nodal layers) 
across the mid-chord boundary layer. 

Figure 6 illustrates some details of the surface 
triangulation, including the “full viscous” and the “wall- 
function” grids. The smooth clustering of surface 
triangulation is evident on the wing and fuselage. Note also 
the spanwise stretching of the grid along the leading edge in 
Fig. 6(b) that reduces the overall number of cells 
considerably. Figures 6(c) and 6(d) illustrate the respective 
extents of the advancing (viscous) layers for the "full 
viscous” and “wall function” grids. 

The provided CAD representation of the DLR-F4 
model has a wing with a blunt trailing edge definition. The 
trailing edge is resolved on the present grids by five grid 
nodes. Wing trailing edge condition (sharp vs. blunt) has 
been known to have influence on location of the upstream 
shock waves and, thus, on the computed drag counts. The 
effect of wing trailing edge on the drag prediction and 
issues regarding the grid resolution at this location are 
among the subjects requiring further in-depth 
investigations. The present grids are restricted to an 
underlying 0-type topology that limits the attainable 
resolution of the trailing-edge wake region. Automatic grid 
adaptation (currently under development) would further 
facilitate the study of grid sensitivity and its effect on 
solution accuracy. 


USM3Dns Flow Solutions 

Flow solutions were computed on the DLR-F4 
configuration with USM3Dns for the conditions shown in 
Table 2 as recommended by the workshop organizing 
committee. All cases are at a Reynolds number of 3, OX 10 
based on the mean aerodynamic chord (mac). Cases 1 and 2 
were designated as “required” computations by all 
participants. The participants were also required to use the 
workshop-provided grids for Case 1 . Ffowever, they were 
encouraged to generate their own computational ^ids or, if 
desired, use the provided grids for the remaining cases. 
Computations for Cases 3 and 4 were optional. An 
additional case was computed with USM3Dns for Mx= 
0.75, a=0.93 deg., and Remac=3.0X 10^ to initially evaluate 


Case 

M. 

a, deg. 

Cl (±0.001) 

1) Single 
point 

0.75 

- 

0.5 

2) Drag 
polar 

0.75 

-3.0, -2.0. 
-1.0, 0.0, 
1.0, 2.0 

- 

3) Mach 
sweep 

0.50, 0.60, 
0.70, 0.75, 
0.76, 0.77, 
0.78, 0.80 

- 

0.5 

4) Drag 
rise 

0.50, 0.60, 
0.70, 0.75, 
0.76, 0.77, 
0.78, 0.80 

- 

0.40, 

0.50, 0.60 


Table 2. Workshop flow cases for the DLR-F4 
configuration at Remac=3.0X10*^ 

the code perfomiance and assess preliminary comparisons 
with Ref 4. All workshop computations assumed fully 
turbulent flow, whereas the experimental data utilized a 
fixed transition technique to trip the flow from laminar to 
turbulent near the leading edge. 

Boundary conditions for the USM3Dns solutions 
consisted of a mirroring condition for the symmetry plane, 
characteristic inflow/outflow for the outer boundaries, a 
wall function on all aerodynamic surfaces with an exception 
of the wing trailing edge. The thick trailing ed^e surface 
was treated by a special “wake” artifice^^ that was 
developed to mimic the relieving effect of a blunt-base 
wake on a coarse grid. This formulation introduces a 
solution-defined transpiration velocity on the blunt-base 
boundary faces to provide a smooth departure of the flow 
past the comer. This boundary condition has been tested 
and is used extensively with inviscid flows and wall 
function applications. 

The USM3Dns flow solutions were converged to at 
least 4-orders of magnitude residual reduction which 
required between 2000 and 5000 iterations, depending on 
the extent of flow separation. All solutions were generated 
with no difficulty. Each flow condition required between 
1.7 and 4.2 wallclock hours to run on 48 processors of the 
NASA Origin 2000 National Aerospace Simulator (NAS) 
computer. 

Figure 7 portrays the surface pressure coefficient 
contours on the wing upper surface for the initial solution. 
At the flow condition 0.75, a=0.93 deg., the surface 
pressures show a classic transonic shock formation with an 
inboard lambda compression region that is consistent with 
the CFD and experimental data reported in Ref 4. 
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Figure 8 depicts a small vortex separation at the trailing 
edge wing/fuselage junction. Ref. 4 notes the experimental 
and computational presence of this feature that is expected 
to have some effect on total drag. It is noteworthy that the 
wall function does not impede the formation of such three- 
dimensional flow separations. 

EfTect of Flux Limiter 

During the analysis of the initial flow solution, the lift 
and drag coefficients appeared high compared to the 
experimental result in Ref 4, as observed in the first two 
rows of Table 3. Other users of the code have also observed 
this tendency with similar configurations in the past. After 
some additional exploratory calculations, the minmod flux 
limiter was isolated as the cause. As evident in Table 3, 
turning off the limiter resulted in a 48-count reduction in 
total drag and a 30-count increase in lift coefficient. 

Flux limiters are routinely used in upwind schemes to 
prevent numerical overshoots in pressure just ahead of 
shock waves. Figure 9 portrays the effect of the limiter on 
predicted surface pressure coefficients (Cp) in comparison 
with experimental data at six span stations on the wing. In 
general the flux limiter has little effect in altering shock 
location, but both tend to predict the shock location slightly 
aft. Note also that the unlimited solution exhibits no sign of 
an overshoot of Cp near the shock, suggesting that the 
limiter is not needed for this application. While the effect of 
the limiter on the Cp distributions appears insignificant, the 
drag differs by 48 counts. Close observation of Fig. 9 
reveals that the limited solution (solid line) has slightly 
more positive Cp at the leading-edge suction peak, and 
slightly more negative Cp on the aft-wing upper surface 
than that for the unlimited solution (dashed line). One can 
infer that a higher pressure acting on the forward- facing 
leading edge surface, and lower pressure acting on the aft- 
facing aft-wing upper surface will each contribute a positive 
axial component to pressure drag. Thus the subsequent 
solutions for the workshop cases have been run without a 
flux limiter due to this unresolved issue. 


Data 

C 

L 

c 

D 

Experiment (Ttef. 4) 

0.602 

0.0352 

USM3Dns 

0.650 

0.0412 

(with limiter) 

USM3Sns 

0.653 

0.0364 

(no limiter) 


Table 3. Lift and drag coefficients for DLR-F4 model. 
Mx= 0.75, a=0.93°, and Rem„=3.0X 10*. 



Figure 7. Surface pressure coefficient contours on 
DLR-F4 wing, 0.75, a =0.93^ 

Remac=3*0X 10^, (Magenta - low pressure. 
Green - high pressure). 



Figure 8. Predicted oil flow traces on the DLR-F4 


configuration showing footprint of a vortex 
separation at wing-body junction. 
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Experimental Data 
USM3Dns (with limiter) 
USM3Dns (without limiter) 



-1.5 


-l.O 


Cp 


-0.5 


0.0 
0.5 
1 0 





Figure 9. Effect of flux limiter on surface pressure distributions on DLR-F4 wing. M^= 0.75, 
0=0.93°, and Re„„=3.0X 10*. 


Statistical Comparison 

A robust statistical evaluation of 35 CFD solutions 
presented in the Drag Prediction Workshop has been 
compiled by Hemsch in Ref. 16. These solutions were 
generated from a range of structured and unstructured 
Navier-Stokes flow codes, with one- and two-equations 
turbulence models, and from different users. Thus, the 
current statistical analysis serves as an initial baseline to 
assess solution accuracy of an individual code relative to a 
larger population. 

The statistical limits to be used in Tables 4 to 7 were 
extracted from Ref. 16. These tolerances generally express 
the statistical probability that an individual solution would 
fall outside these limits more than one time in 100 in the 
long run by chance (99% coverage). 

The first workshop flow Case 1 is a single point 
computation at - 0.75, Cl ^ 0.5, and Rcmac = 3.0 X 10^. 
Table 4 displays the lift, drag, and pitching moment 
coefficients for two USM3Dns results and from three 


different wind tunnels. A statistical composite of the wind- 
tunnel data is included from Ref. 16 to convey the scatter in 
the data. The electronic data tables that provided the wind 
tunnel coefficients were furnished to an accuracy of 10 
counts. For the purposes of the workshop, the round-off 
error was reduced considerably by use of an interpolation 
technique developed by Vassberg". The USM3Dns flow 
computations were performed at a constant Cl of 0.5 using 
an automated feedback loop in the code that adjusted angle 
of attack to -0.300 deg. In addition, solutions were obtained 
at an averaged angle of attack corresponding to Cl=0.5 
from the three wind tunnel tests (aav^=0.175 deg.). The 
difference of the USM3Dns point result in Table 4 from the 
mean value of the workshop CFD solutions is enclosed 
within parenthesis. While the variation for the computed 
drag coefficient at “constant Cl” (-0.0009) is only 1 count 
outside of the experimental scatter (±0.0008), the angle of 
attack and pitching moment coefficient are significantly 
farther outside their respective ranges. 

More understanding of these differences can be derived 
from the Cp comparisons in Figure 10. Shown is a 
comparison of the two computational solutions in Table 4 
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with the three sets of experimental data at six span stations. 
When Cl is matched between USM3Dns and experiment by 
varying angle of attack (solid line), the distributions differ 
considerably, while the integrated lift is the same. 
Similarly, if an averaged angle of attack is used (dashed 
line), the general agreement with experiment is better ahead 
of the shock, but the shock position is too far aft, and lift 
and drag are higher. Similar mismatch of the predicted 
surface pressures was experienced by most of the other 
workshop participants as well. The many reasons offered at 
the workshop for these differences concerned offsets in 
angle of attack, wind tunnel walls and blockage effects, 
sting effects, aeroelastic effects, and the influence of the 
wing trailing-edge bluntness on the upstream flow. More 
discussion on this matter is available in Ref. 2. As a general 
rule, these comparisons suggest that the experimental data 
set should only be considered as a guide in this study due to 
the number of issues clouding the correlations. The more 
useful focus of our assessment should be on comparisons to 
the statistical composite of CFD solutions from the Drag 
Prediction Workshop. 

A comparison of predicted drag from USM3Dns to the 
larger population of CFD codes is presented in Table 5, and 
similarly in Table 6 for angle of attack and pitching 
moment. The total drag estimate in Table 5 is within 17- 
percent of the scatter of the larger population of codes. 
Angle of attack in Table 6 is within 12-percent and pitching 
moment within 9-percent. These results confirm that the 
accuracy for predicting drag and other aerodynamic 
parameters at transonic cruise conditions using a tetrahedral 
cell-centered scheme coupled with a wall function is well 
within the scatter band of other structured and unstructured 
codes solving the Navier-Stokes equations to the wall. 

Case 2 designates a drag polar prediction at M«, = 0.75 
and six angles of attack: -3, -2, -1, 0, I, and 2 degrees. 
Plots of angle of attack vs. lift, drag vs. lift, and lift vs. 
pitching moment coefficients are presented in Figure 11. 
There is a good consistency among the lift and drag 
measurements in three different wind tunnels. While the 
computational Cp vs. Cl curv^e matches the experimental 
data fairly well, the lift coefficients are generally over- 
predicted at individual angles of attack. The slope of the a 
vs. Cl curve, however, matches those of the experiments 
reasonably well before it breaks downward prematurely at 
higher incidence. As illustrated in Fig. 11, the entire 
predicted a vs. Cl curve appears to lag the experiments by 
an incidence offset of about one-half degree. Similar over 
prediction of the lift has been observ ed from the majority of 
workshop computations. As stated in Ref 2, an average 
offset of 0.414 degree was required to achieve the target Cl 
by all computations. 

Unlike the force measurements, the experimental 
pitching moments in Fig. 1 1 exhibit a larger disparity (at 
the scale presented), and the computational results are 
generally under-predicted (nose-down). While the exact 
source of inaccuracy is not known, several factors are 
believed to have contributed to the mismatch. Among those 
suggested in Ref 2 are the free-air assumption for 
computations, absence of the sting mount in the CFD 


Data 

a 

(deg.) 

CL 

CD 

(Interpolated) 

CM 

ONERA 

0.192 

0.50 

0.0289 

-0.126 

NLR 

0.153 

0.50 

0.0288 

-0.130 

DRA 

0.181 

0.50 

0.0281 

-0.137 

Statistical 

0.18 

0.50 

0.0286 

-0.130 

Composite 

±0.04 

±0.0008 

±0.01 

USM3Dns 

-0.300 

0.50 

0.0277 

-0.158 

(constant Cl) 

(-0.48) 

(-0.0009) 

(-0.028) 

USM3Dns 

0.175 

0.56 

0.0303 

-0.156 

(constant a) 

f-0.005) 

(+0.06) 

(0.0017) 

(-0.026) 


Table 4. Forces and moment data for Case I DLR-F4, 



Cd 

pressure 

Cd 

Skin 

Friction 

Cd 

Total 

Statistical 

Composite 

Experiment 


- 

0.0286 

±0.0008 

Workshop 

Casel 

0.0166 

±0.0036 

0.0134 

±0.0039 


Point 

Result 

USM3Dns 

Casel 


0.0118 

(-0.0016) 



Table 5. Comparison of USM3Dns drag prediction to 
statistical data'^. 



a, deg. 

Cm 

S3 .'S 

U VI 

o 

Experiment 

0.18±0.04 

-0.13±0.01 

o 

cj 

Workshop 

Casel 

-0.26±0.34 

-0.160±0.02I7 

Point 

Result 

USM3Dns 

Casel 

-0.300 

(-0.04) 

1 -0.1581 

(0.0019) 


Table 6. Comparison of USM3Dns angle of attack 
and pitching moment to statistical data^^. 
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Figure 10. Surface pressure distributions on DLR-F4 wing. 
Case 1: M.. = 0.75, Cl = 0.5, and Re„,„ = 3.0 X 10*, 
and M. = 0.75, a = 0.175°, and Re„„ = 3.0 X 10*. 


model, aeroelastic effecls under load, and the absence of a 
tail that affects stability of the model. Although the 
experimental data have been corrected to free-air and the 
CFD model adjusted for aeroelastic effects, these 
corrections are believed to be crude and may have been 
inadequate to warrant similar conditions for the 
experiments and computations. The nose-up break in the 
pitching moment is slightly premature as noted for the lift 
coefficient. This result is unfortunate, but not surprising 
because of the many factors cited in Ref 2. 

Several polar related slope parameters, such as lift- 
curve slope (Ct,a), lift at zero angle of attack (Cl,o). drag at 
zero lift (Cd@Cl= 0), and span efficiency factor, 
( k=(CirCD@CL=0)/(CL)^ ) are compared against workshop 
CFD results in Table 7. The difference of the USM3Dns 
point result in Table 7 from the mean value of the workshop 
CFD solutions is enclosed within parenthesis. The 
USM3Dns results are well within the scatter for each of 



CL,a 

Cl,o 

Cd@ 

Cl=0 

k 

Statistical 

Composite 

Expt 

0.114 

±0.0038 

0.473 

±0.006 

0.0188 

±0.0010 

0.0374 

±0.0080 

Wkshp 

Case2 

0.120 

±0.0067 

0.531 

±0.054 

0.0198 

±0.0041 

0.0365 

±0.0046 

Point 

Result 

USMlDns 

Case2 

0.123 

(+0.003) 

0.539 

(+0.008) 

0.0190 

(-0.0008) 

0.0337 

(-0.0028) 


Table 7. Comparison of USM3Dns lift-cur\e slojge and 
drag-polar parameters to statistical data”^. 


11 


American Institute of Aeronautics and Astronautics 







o 

ONERA Data 

□ 

DRA Data 

o 

NLR Data 


USM3Dns 


Figure II, Force and moment data for DLR-F4, Case 2: 
M„ = 0.75, and Re„,.c = 3.0 X 10‘. 


these parameters, and generally within the scatter of the 
experimental composite. This suggests that CFD data can 
be as useful as wind tunnel data for studying trends in 
aerodynamic forces. 

The last two optional Cases 3 and 4 are representative 
of constant Cl drag rise predictions for eight different Mach 
numbers and three values of Cl- The Mach sweep, drag 
rise curves are useful for analyzing the drag characteristics 
of commercial aircraft and are often used by the industry. 
Figure 12 illustrates the predicted drag rise curves, along 
with the experimental data, for Cl values of 0.4, 0.5, and 
0.6. At the Cl of 0,4, the computation agrees with the 
experiments reasonably well, however, the agreement 
degrades with increasing lift and Mach number where 
additional aeroelastic defonuation and separation onset 
occur. Recall that the computational grid was constructed 
on a single-point surface estimate of aeroelastic 
deformation. Additional load-induced grid deformations are 
needed for more accurate prediction of aerodynamic 
coefficients across the Mach range. The prediction of 
separation onset is generally difficult due to sensitivities 
associated with transition and turbulence models. Although 
the drag is generally under predicted at these conditions, the 
onset of the drag rise has been properly captured by the 
current computations. As in the preceding cases, the 
mismatch in drag rise between CFD and experiments was 

common among the workshop CFD results. 


CONCLUDING REMARKS 


The specifications of three unstructured tetrahedral 
grids that were generated for the first ATAA Drag 
Prediction Workshop held at the 19^^ ATAA Applied 
Aerodynamics Conference in Anaheim, CA, 9-10 June 
2001 are presented. The underlying grid generation 
methodology embodied in the VGRIDns code is also 
reviewed. The VGRIDns code is part of a complete flow 
analysis system known as TetrUSS. Special emphasis is 
placed on the utility of this system as reflected by ease of 
use, and total turnaround time from CAD to final 
postprocessed drag polar solutions in days. 

An assessment of the TetrUSS flow solver USM3Dns 
for predicting drag on a transonic transport configuration is 
made relative to three experimental data sets and a 
statistical composite of 35 CFD results from a range of 
structured and unstructured flow solvers using one- and 
two-equation turbulence models. Over 30 USM3Dns flow 
solutions on the DLR-F4 configurations were generated for 
this study with no experienced difficulty. Each solution 
required between 1.7 and 4.2 wallclock hours to run on 48 
processors of the NASA Origin 2000 NAS computer. The 
assessment led to the following conclusions: 
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1 . Initial studies exposed a deficiency in the USM3Dns 
flux limiter and demonstrated that the limiter may not 
be needed for the present class of applications. 

2. The available DTR-F4 experimental data set should 
only be considered as a guide in this study due to the 
number of issues clouding the correlations. The more 
useful focus of an assessment should be on 
comparisons to the statistical composite of CFD 
solutions from the Drag Prediction Workshop. 

3. A comparison of predicted drag from USM3Dns to the 
larger population of CFD codes has demonstrated that 
its accuracy for predicting drag at transonic cruise 
conditions using a tetrahedral cell-centered scheme 
coupled with a wall function boundary condition is 
well within that of other structured and unstructured 
codes solving the Navier-Slokes equations to the wall. 
The wall function is beneficial for reducing solution 
stiffness and the required number of cells across the 
boundary layer. 

4. Several polar related slope parameters, such as lift- 
curve slope, lif\ at zero angle of attack, drag at zero lift, 
and span efficiency factor are compared against 



workshop experimental and CFD results. USM3Dns 
results are w'ell within the scatter for each of these 
parameters, and generally within the scatter of the 
experimental composite. This suggests that CFD data 
can be as useful as wind tunnel data for studying trends 
in aerodynamic forces. 

5. The combination of rapid grid generation and flow 
solution turnaround exhibited by VGRIDns and 
USM3Dns, respectively, along with the demonstrated 
accuracy for predicting aerodynamic forces and 
moments as compared to other structured and 
unstructured solvers underscores the practicality and 
usability of the TetrUSS system. 

The comparison of USM3Dns solutions to a statistical 
composite from a larger body of CFD codes has been very 
beneficial. As a result, several key issues were surfaced 
that merit further consideration in future workshops. 
Among these are static aeroelasticily, use of flux limiters, 
flow transition, and grid resolution in the trailing-edge 
wake region. 



MACH 


o 

ONERA Data 

□ 

DRA Data 

o 

NLR Data 


USM3Dns 


Figure 12. Drag rise curves for DLR-F4. 
Cases 3&4: Cl= 0.4, 0.5, and 0.6, 
RCmac = 3.0 X 10*. 
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